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Sophistication of microwave radiative transfer models 



= /(radiative transfer interactions in the snowpack) 

RT models relate snow and radiation properties 


Models can be: 


Theoretical 


Semi empirical 


Empirica 
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D M RT- M L - Snowpack and Soil Configurations 

(Dense Media Radiative Transfer Multi-Layer model) 
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Every layer is defined by: 

- Grain size (sphere radius r) 

- Density (p) 

- Temperature (T) 

- Stickiness parameter (r) 

- Liquid water content (LWC) 
It can be layers of snow or 
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D M RT- M L - Snowpack and Soil Configurations 

(Dense Media Radiative Transfer Multi-Layer model) 
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It can be layers of snow or 


Underneath interface can be: 

Snow, Ice, Soil (smooth or rough). Fresh water 
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D M RT- M L - Snowpack and Soil Configurations 

(Dense Media Radiative Transfer Multi-Layer model) 


Interface model 
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Dielectric constant 

ID 
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SM: soil moisture (volume fraction) 

e: dielectric constant 

cr: surface RMS height (m) 

/clay- 3^d /sand- fractions of clay and sand 
Porga: density of dry organic matter (kg m“^) 
T: soil temperature (K) 

Q, and H: dimensionless parameters 


P99 = Pulliainen et al. (1999) 
D85 = Dobson et al. (1985) 
M87 = Matzler (1987) 
M06 =Matzler et al. (2006) 


p. 6 ghislain.picard@ujf-grenoble.fr & ludovic.brucker@nasa.gov The DMRT-ML model 


Context DMRT-ML Validations Challenges Conclusion 

D M RT- M L - Snowpack and Soil Configurations 

(Dense Media Radiative Transfer Multi-Layer model) 
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Every snow layer is defined by: 

- Grain size (sphere radius r) 

- Density (p) 

- Temperature (T) 

- Stickiness parameter (r) 

- Liquid water content {LWC) 


Underneath interface can be: 

Snow, Ice, Soil (smooth or rough). Fresh water 
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DMRT-ML - Snowpack and Soil Configurations 

(Dense Media Radiative Transfer Multi-Layer model) 
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Every snow layer is defined by: 
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Underneath interface can be: 

Snow, Ice, Soil (smooth or rough). Fresh water 
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D M RT- M L - Elements 
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DMRT-ML computes: 

. emission 'j 

. scattering > using the DMRT theory 
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DMRT-ML computes the reflection coefficients at: 
. snow/snow interfaces 
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using Fresnel 


DMRT-ML solves the radiative transfer equation 

using the DISORT method 
(64 or 128 streams are adequate) 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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D M RT- M L - Elements 


Optional correction for large particles (modified from Grody, 2008)\ 


the scattering efficiency coefficient cannot exceed 2 

^Qmax_2^ theoretical asymptotic value for very large particles) 



Frequency (GHz) 


It is not recommended to consider a many and thick layers of large grains _ 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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D M RT- M L - Elements 
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. QCA-CP 

. Rayleigh assumption 
. No large particles 

. Poly-disperse 
(i.e. Rayleigh distribution) 
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DMRT-ML - Implementation 
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D M RT- M L - Open source model 


DMRTML model (Multi-layered Radiative Transfer Model For Passive Microwave) / DMRT-ML 


DMRTML model (MuttL.. X 


® lgge. 0 BugJr/~picard/dmrtml/ 


It dmrt-mi 


a 


Obtaining the code: 


version 1.6 


The source code and the technical documentation can be downloaded at 
ihttp ://lgge.osug.fr/ ~pi c ardAimilinl/do a dA 

Please read and accept the license: GNU General Public License version 3 
(http://www.opensource.Org/licenses/gpl-3.0.html). In particular, if you modify the 
code, it must remain open source with the same GPL-3.0 bcense. Also note that 
this software comes with no warranty. 


Mailmg list and (optional) registration: 


To get help and receive updates about DMRT-ML, or discuss about appbcations of 
passive microwave remote sensing, please subscribe to the obs-dmrtniL 
mailing-list . 

In addition, it is recommended but optional to register. It helps us to keep track 
of DMRT-ML users and provides justification for funding! ff you want your 
publications using DMRT-ML to appear on this page, use this form or post a 
message on the mailing-bst. 


FiiU name: 

[ ^^1 

Email: 

f 1 

Institution 

and 

country: 




Your 

research 

activity: 



1 RegBter | 


Installation and compilation: 
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D M RT- M L - Open source model 
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► Class List 
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m od_d m rtparamete rs 


med_soil 
Module Members 
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F ur>ctk>nrs/SubroutiTies 


mod_dmrtparameters Module Reference 

Compute dense media radiative transfer parameters [extinction and single scattering albedo) under various assumptions [monodispese or polydisperse, Grody’s limitation^ More.. 

Functions/Subroutines 


subroutine dmrtparameters_grcH}yapproach [trequency, temperature^La.tau, LWC, medium^ EefflD, Eeff, albedo, beta, Eice) 

Compute radiative transfer parametersjnrti^tiDn and single s(»itt|ring albedi) using DWRT theory under the approximations: QCA-CR small pagicules, sticky short range, single species. In 
addition, an empihcaJ correction fer larle p^cu^;Japi;ji^ tmsld acLGufm |i rnyd Rla'VA | ^ E-^i [rngi h ■ it iUs still far from perfect tts main 

purpose is to prevent DMRT divergenclfor lwg»kilul^ihln aye r^o ilagicy rli n i hnvfv^, jiM L|di|o tie iSlidlrerjtipcJ to investigate the scattering 

behavior of large particules. To disablelhii^-lVj^ABffclule li^MfaiVrtwe dmwpAaWrerViastlsAfriiJQni.^^ I I w wJ w I I I 


subroutine dmrtparametera [frequency, temperature, f, a, tau, LWC, medium, EeffO, Eeff, albedo, beta. Bee) 

Compute radiative transfer parameters [extinction and single scattering albedo) using DMRT theory under the assumption: QCA-CR, small particules, sticky short range, single species. 

subroutine dmrtparameters_dlst [frequency,temperature,f,a,tau, LWC, medium, EeffO, Eeff, albedo, beta, Eice) 

Compute radiative transfer parameters [extinction and single scattering aJbedo) using DMRT theory underthe assumption s: QCA-CR small particules,* no ne-sticky*, multispecies. Grain size 
repartition is given by a Rayleigh distribution. 


Detailed Description 

Compute dense media radiative transfer parameters [extinction and single scattering albedo) under various assumptions [monodispese or polydisperse.Grody’s limitation. 


Function/Subroutine Documentation 


subroutine mod_ dm rtparamete rsiidmrtparameters { real*3,intent[in) 

reai*Bpintent[in) 
real*&pintent(in) 
real*Bpintent[in) 
rears Jntent[in) 


frequency r 
temperature, 
f, 
a, 

tau, 


Gefwmted on Wed Apr 11 2012 08:42:43 for DMRTML by ffllO 


The one-line Python code: 

dmrtml.dmrtm I (frequency, 128, height, density, radius, temp, 
tau=dmrtml.NONSTICKY,dist=False,soil^=None) 


■^c\(y 
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Methodology 


Measurements 

T(z,t), p(z) & r^p,(z) 


DMRT-ML 


Brightness Temperature 

"*”6 ig(h & "*"537(1) 


ropt may be measured using near infrared sensors 
(camera, integrating sphere, profiler, etc.) 


□ I ^ ■^c\(y 
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Methodology 


Measurements 

T(z,t), p(z) & r^p,(z) 


1 


DMRT-ML 



Brightness Temperature 

“► Tg ^g(t) & Tgg^(t) 



► RMSE 


Tg AMSR-E 2007 
19 and 37 GHz 


q; is a calibrated parameter identical at both frequencies. 


(Brucker et a!., 2011 JoG) 
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Dome C, Antarctica 


Measurements 

T(z,t), p(z)&r^pj(z) 


DMRT-ML 


Brightness Temperature 

— ► 19 W ^ 3y(^) 





RMSE 


Calibration result: a= 2.8 


Tg AMSR-E 2007 
19 and 37 GHz 



V pol. 

RMSEig ~0.4K 
RMSE 37 ~ 1 K 


RMSE -0.8K 


(Brucker et a!., 2011) 
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Dome C, Antarctica 


Simulated and observed Tb at 37 GHz, H pol. 



Rapid variations are not reproduced (constant surface properties) 


(Brucker et a!., 2011 JoG) 
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Dome C, Antarctica 


Aquarius (1.4 GHz) TB observations 
show rapid variations at H pol. 



Date 


29 . 2 ° 

38 . 4 ° 

46 . 3 ° 
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Dome C, Antarctica 


Measurements 

T(z), p(z) & rjz) 


DMRT-ML 



Brightness Temperature 


Tbi9 & Tg 3 ^ 


Tg SBR 


19 and 37 GHz 


-► RMSE 


Calibration result: cx=2.3 



Reasonable simulations at both 
Vertical Si Horizontal pol. 

Reproduce smooth and chaotic 
surfaces 

Issues at high incidence angles 
(> 60°) 


(Picard et a!., 2014 TC) 
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Ablation Zone, Antarctica 


Cap Prud’Homme 

(Dupont et al., 2014 TGRS) 
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Ablation Zone, Antarctica 


Cap Prud’Homme 

(Dupont et al., 2014 TGRS) 



Bubbly Ice 



Importance of 
bubble size S>l p-,ce 


The DMRT theory is 
valid up to fractional 
volume of 30%. 
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Ablation Zone, Antarctica 


Cap Prud’Homme 

(Dupont et al., 2014 TGRS) 



Bubbly Ice 



Importance of 
bubble size S>l p-,ce 


The DMRT theory is 
valid up to fractional 
volume of 30%. 
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Ablation Zone, Antarctica 


Cap Prud’Homme 

(Dupont et al., 2014 TGRS) 



Bubbly Ice 



Importance of 
bubble size S>l p-,ce 


The DMRT theory is 
valid up to fractional 
volume of 30%. 
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Ablation Zone, Antarctica 


Cap Prud’Homme 

(Dupont et al., 2014 TGRS) 
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37 H-pol (K) Simulated (K) 


Context DMRT-ML Validations Challenges Conclusion 

Seasonal Arctic/sub-Arctic Snowpacks: Manitoba Quebec 
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Context DMRT-ML Validations Challenges Conclusion 

Dependence between H V pol. 


Observed and simulated TB 

over Dome C (vr), ice-sheet ablation or percolation areas (o), 
tundra (A), windy tundra (v)- fen (D)- grassland (o) 



DMRT-ML predicts reliable dependence between H V pol. near 50-55^ 
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Context DMRT-ML Validations Challenges Conclusion 

Validations Si Applications 


DMRT-ML has been validated, in very different environments, on 

- deep firn on ice sheet 

- shallow snow covers overlying soil (Arctic and Alpine regions) 

or ice 
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Context DMRT-ML Validations Challenges Conclusion 

Validations Si Applications 


DMRT-ML has been validated, in very different environments, on 

- deep firn on ice sheet 

- shallow snow covers overlying soil (Arctic and Alpine regions) 

or ice 


DMRT-ML is now used in applications for: 

- hemispheric SWE retrievals 

- sensor developments 

- snow property retrievals with inverse methods 

- coupling with snow evolution models 

- data assimilation 
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Context DMRT-ML Validations Challenges Conclusion 

Outline 


1. Descriptions of the DMRT-ML Model 

2. Validations Experiments 

3. Current Limitations Challenges 

4. Conclusion 
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Current Limitations ^ Challenges 


Liquid water content 



Water is in the top 10 cm 
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Current Limitations Challenges 

Grain Size Factor (a) 


Density 


Roughness 


Large Particles 


□ ; ^ -El -E -E '00^0 
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Current Limitations ^ Challenges 

Grain Size Factor (a) 


. a, factor between r^l^surec/ ^ is likely due to microstructure. 

. Stickiness and grain size distribution may contribute to a. 

Density 


Model 

Location 

a 


Antarctica 

2.8 

DMRT-ML 

Antarctica 

2.3 

Barnes Ice Cap 

3.5 


Sub-Arctic 

3.3 


Roughness 


Large Particles 
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□ I 

The DMRT-ML model 







Scattering and absorption coefficients (nrr‘) 
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Current Limitations Challenges 


Grain Size Factor (a) 

Model 

Location 

a 



Antarctica 

2.8 



DMRT-ML 

Antarctica 

2.3 



Barnes Ice Cap 

3.5 




Sub-Arctic 

3.3 



. a, factor between r^l^surec/ ^ is likely due to microstructure. 

. Stickiness and grain size distribution may contribute to a. 


Density 

. The DMRT theory is valid up to fractional volume of 30%. 


i.e. 0 - 275 kg m 
and 642 - 917 kgm“^. 

Bridging = Combination of the two 

valid domains (Dierking et al., 2012) 

''"'0 100 200 300 400 500 600 700 800 900 1000 □ 

Densltyjkg 
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Current Limitations Challenges 

Grain Size Factor (a) 

. a, factor between Su. is likely due to microstructure. 

. Stickiness and grain size distribution may contribute to a. 

Density 

. The DMRT theory is valid up to fractional volume of 30%. 

Roughness 

. Surface roughness, and layer interface roughness are not considered. 

Large Particles 


Model 

Location 

a 


Antarctica 

2.8 

DMRT-ML 

Antarctica 

2.3 

Barnes Ice Cap 

3.5 


Sub-Arctic 

3.3 
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Current Limitations Challenges 

Grain Size Factor (a) 


. a, factor between r^|3sured ^ is likely due to microstructure. 

. Stickiness and grain size distribution may contribute to a. 

Density 

. The DMRT theory is valid up to fractional volume of 30%. 

Roughness 

. Surface roughness, and layer interface roughness are not considered. 

Large Particles 

. The DMRT-Mie theory is only available under the QCA 

is time consuming 

□ I ^ 
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Model 

Location 

a 


Antarctica 

2.8 

DMRT-ML 

Antarctica 

2.3 

Barnes Ice Cap 

3.5 


Sub-Arctic 

3.3 
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Current Limitations ^ Challenges 

Brine in snow and ice 


DMRT-ML does not include 
dielectric constant parameterizations for 
brine content (yet) 


Configuration already possible with DMRT-ML 
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Context DMRT-ML Validations Challenges Conclusion 

Outline 


1. Descriptions of the DMRT-ML Model 

2. Validations Experiments 

3. Current Limitations Challenges 

4. Conclusion 
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◄ 4 ^ k < -E > -E '0<\(y 


The DMRT-ML model 
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Conclusion 

DMRT-ML is validated and operational for both seasonal snow 

& perennial snow 


DMRT-ML is a published open source model available at: 

http: / /Igge. osug.fr/ '^picard /dmrtml / 

Email distribution list obs-dmrtml@ujf-grenoble.fr 


DMRT-ML Community Model 


Contributions for model developments are welcome 
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- A. Roy, G. Picard, A. Royer, B. Montpetit, F. Dupont, A. Langlois, C. Derksen and N. Champollion, Brightness 
temperature simulations of the Canadian seasonal snowpack driven by measurements of the snow specific surface 
area, IEEE TGRS, vol.51, no. 9, pp. 4692-4704, 2013. 

- F. Dupont, G. Picard, A. Royer, M. Fily, A. Langlois, A. Roy, and N. Champollion, Modeling the microwave 
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vol.52, no. 10, pp. 6639-6651, 2014 

- G. Picard, A. Royer, L. Arnaud, and M. Fily, Influence of snow dunes on the microwave brightness temperature: 
investigation with ground-based radiometers at Dome C in Antarctica, The Cryosphere, vol.8, 
doi:10.5194/tc-8-l-2014 

Applications 


- L. Brucker, G. Picard, and M. Fily. Snow grain size profiles deduced from microwave snow emissivities in 
Antarctica. J. of Glaciol., 56(197), 2010. 
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